Bulk Aluminum at High Pressures: A First-Principles Study 
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The behavior of metals at high pressure is of great importance to the fields of shock physics, 
geophysics, astrophysics, and nuclear materials. We study here bulk crystalline aluminum from first- 
principles at pressures up to 2500 GPa - soon within reach of laser-based experimental facilities. Our 
simulations use density-functional theory and density-functional perturbation theory in the local- 
density and generalized-gradient approximations. Notably, the two different exchange-correlation 
functionals predict very similar results for the fee — > hep, fee —> bee, and hep — > bee transition 
pressures, around 175 GPa, 275 GPa, and 380 GPa respectively. In addition, our results indicate that 
core overlaps become noticeable only beyond pressures of 1200 GPa. From the phonon dispersions 
of the fee phase at increasing pressure, we predict a softening of the lowest transverse acoustic 
vibrational mode along the [110] direction, which corresponds to a Born instability of the fee phase 
around 725 GPa. 



First-principles calculations have proved useful to the 
fields of geophysics, 1 astrophysics^ and nuclear materi- 
als^ Aluminum, being cubic close-packed and having no 
d-shell electrons, is a prototype for theoretical predic- 
tions and understanding the high-pressure behavior of 
simple metals^ Currently the National Ignition Facility 5 
at LLNL is expected to achieve shockless compression 6 
of metals up to 2000 GPa. This new facility may provide 
rapid advancements to high-pressure physics and could 
partner very successfully with theoretical studies. 

The equation of state (EOS) and phase stability of 
aluminum were first studied from first-principles in the 
early 1980sji&2 In all cases the predicted phase se- 
quence was fee — > hep — > bee, but predictions differed 
greatly in the transition pressures. Several other calcula- 
tions within the local-density approximation (LDA) 1 ^ or 
the generalized-gradient approximation (GGA) 11 ' 12 have 
since then been performed, with a predicted static (i.e. 
without the phonon contribution) fee — > hep transition 
pressure of 205 ± 20 GPa 10 in LDA and 170 GPaii and 
192 GPai2 in GGA. These discrepancies are more notable 
for the hep — > bee transition pressure: 565 ± 60 GPai°- in 
LDA versus 360 GPa 11 in GGA, leaving significant uncer- 
tainties open. Theoretical work on the vibrational prop- 
erties of aluminum also suggests for the fee — > hep transi- 
tion a transition pressure higher than the static one. 11 ' 12 
Elastic propertie d 13 ' 14 and the absolute strength under 
tension 1 ^ have also been calculated; the latter results are 
of particular interest as they demonstrate the important 
role vibrational modes play in determining mechanical 
stability and suggest that shear failure modes are inher- 
ent in aluminum. 

Experimentally, the equation of state of aluminum at 
high pressures was studied by shock-compression^ 6 - at 
pressures above the predicted maximum for the fee — > 
hep phase boundary^ but a transition was not observed. 
However recent diamond anvil cell experiments observed 
a fee — » hep transition at 217 ± 10 GPa 17 highlighting 
the difficulty in achieving thermodynamic equilibrium in 
shock-compression. 

In this article we report first-principles calculations 



of aluminum under hydrostatic compression up to 2500 
GPa. In order to assess mechanical stability under shock, 
we also calculate the vibrational properties in the fee 
phase and determine the elastic constants from the slopes 
of the phonon dispersions (i.e. the sound velocities). 

The equations of state in the fee, bec, and hep 
phases have been calculated with density-functional the- 
ory (DFT) within both LDA±& and GGA^ Calculations 
have used the Quantum-ESPRESSO package. 20 We 
use plane-wave basis sets and pseudopotentials and both 
3 electron (3e) norm-conserving pseudopotentials,— with 
the 3s and 3p electrons in the valence and nonlinear core- 
corrections, and 11 electron (lie) ultrasoft pseudopoten- 
tials^ 2 - where the 2s and 2p electrons, usually frozen in 
the core, are included explicitly in the valence. The in- 
clusion of the 2s and 2p electrons in the valence is es- 
sential to investigate the relevance of inner core electrons 
at high pressure. The planewave cutoffs for the wave- 
functions are 25 Ry and 100 Ry for the 3e and lie pseu- 
dopotentials, respectively, and 150 Ry and 800 Ry for 
the charge density. Brillouin zone integrations have been 
performed using a cold smearing 14 of 0.02 Ry over shifted 
Monkhorst-Pack meshes of order 16x16x16 for the fee, 
22 x 22 x 22 for the bee and 16 x 16 x 10 for the hep 
phases. The large sizes of the k-point meshes are neces- 
sary to obtain fully-converged transition pressures. The 
data for the total energy as function of volume have been 
fitted to the Birch 3rd order EOS^ 3 - near equilibrium to 
obtain equilibrium volumes and bulk moduli. The data 
for the equation of state have been determined from cal- 
culations performed at around 50—100 different volumes 
over a pressure range of — 2500 GPa. Finally, the vi- 
brational properties have been computed using density- 
functional perturbation theory (DFPT).— The dynami- 
cal matrices have been calculated on a 4 x 4 x 4 g-point 
mesh and Fourier interpolation has been used to evaluate 
the phonon frequencies on finer grids. 

The role of the inner core electrons is of primary con- 
cern at very high pressures. Under normal conditions, 
there is not sufficient overlap between the core and va- 
lence shell electrons to question the frozen-core approxi- 
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TABLE I: The equilibrium lattice parameters, bulk moduli, 
and transition pressures of Al calculated with the different 
pseudopotentials described in the text, and compared to ex- 
perimental results. 



(a) 



a [A] Bo [GPa] 


fee — > hep fee — » bee hep — > bee 


lie GGA 4.044 73.2 


175 GPa 275 GPa 383 GPa 


lie LDA 3.985 83.7 


172 GPa 272 GPa 380 GPa 


3e GGA 4.055 71.8 


180 GPa 285 GPa 420 GPa 


Expt. 4.025 a 78.3 a 


217 GPa" - 



(a) 



„ o 







E 

* — i 


-20 


o 




s 


-40 


X 




i 

X 


-60 


ii 

X 


-80 


<3 






-100 



' ' ' ' 








FCC(3e/11e) 






3e BCC 






--- lie BCC 


j i 





500 1000 1500 2000 
Pressure [GPa] 



2500 



(b) 



X! 


12 - 


>, 




E 


8 - 




4 - 


X 




X 


- 


X 


-4 ■ 


<1 


-8 ■ 



— I 

FCC ' 

- - HCP 

BCC 



-L. 



100 200 300 400 
Pressure [GPa] 



500 



FIG. 1: (a) Enthalpy of the bec phase (3e and lie GGA) 
relative to the enthalpy of the fee phase (3e and lie GGA, 
respectively), (b) Enthalpies for hep and bec phase relative 
to the enthalpy of the fee phase. The fee —* hep, fee — > bec, 
and hep — > bee transition pressures are 175 GPa, 275 GPa, 
and 380 GPa respectively. 



mation^ but at the pressures considered here core over- 
laps may become significant. To study the validity of the 
frozen-core approximation, we first compared the equa- 
tions of state for different phases using both the 3e and 
lie pseudopotentials. We report in Table U the equilib- 
rium lattice parameters and bulk moduli in the fee phase 
at zero pressure, and in Fig. [T^i the relative enthalpies 
of fee and bec Al with respect to the fee phase, up to 
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FIG. 2: (a) The c/a ratio of the hep phase plotted with in- 
creasing pressure, (b) Volume versus pressure for the fee, hep, 
and bec phases plotted relative to the fee phase. 



2500 GPa. These results show that 1) there is little dif- 
ference between the LDA and GGA predictions, hinting 
at a broad applicability of density-functional theory in 
either approximation, and 2) that the role of the core 
electrons starts to become noticeable only around 1200 
GPa, even if already at zero pressure the cores of the 
3e pseudopotential start to overlap.— The equations of 
state for aluminum in the fee, bec and hep phase, using 
the lie GGA pseudopotentials, are shown in Fig.[T|D. Al- 
though the lie and 3e calculations give consistent results 
up to 1200 GPa, the calculated transition pressures can 
vary, particularly for the hep — > bee transition. This 
could easily derive from the fact that enthalpy differ- 
ences between these three phases are only a few mRy 
(see Fig. [lb), and so, even at full computational conver- 
gence of all parameters, small effects (e.g. core-state re- 
laxations), which could shift the calculated enthalpy by 
less than a mRy, can significantly affect the calculated 
transition pressures. On the other hand, core electrons 
seem to have a negligible effect in determining the equilib- 
rium volume, bulk modulus, and even phonon dispersions 
(see below). 

Our LDA and GGA results are consistent with previ- 
ously reported GGA calculations ; 11 ^ 2 discrepancies arise 
with the LDA results reported in Ref. [l(| that predict 
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FIG. 3: The frequencies of the two non-degenerate acoustic 
phonon modes at X calculated with both the 3e and lie pseu- 
dopotentials and plotted as a function of increasing pressure. 



205 ± 20 GPa and 565 ± 60 GPa for fee -> hep and 
hep — + bee transition pressures, respectively. This dis- 
crepancy could arise from Ref. [T3 using only 10 — 15 
points to fit the equation of state: as reported there, 
this approximation could significantly affect transition 
pressures due to the aforementioned small enthalpy dif- 
ferences between competing structures. We also observe 
that all parameters of the calculation, and particularly 
the k-point sampling of the Brillouin zone, need to be 
carefully converged. 

Although the pressure that we obtain for the fee — > 
hep transition, 175 GPa, is consistent with previous 
work a 11 ' 12 , this result is lower than the experimental 
value of 217 GPa. As suggested in Refs. [n] and[l2l this 
discrepancy could arise from excluding the phonon con- 
tribution to the free energy - a hypothesis that should be 
thoroughly tested, but that is beyond the scope of this 
brief report. 

It should be noted that in our simulations the hep 
phase was always fully relaxed to identify the optimal, 
equilibrium c/a ratio; this is shown in Fig.[2ji. Compar- 
ison with experiment at 292 GPa finds agreement in the 
c/a ratio to within 0.1% and well within experimental 
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FIG. 4: (a) The lowest energy branch of phonon dispersion 
for fee Al with increasing pressure, (b) As above, enlarged 
around F, in the K — > F direction. A steady flattening is 
observed with increasing pressure. 



five random distinct perturbations. The bec structure 
was always found. 

In order to estimate the dynamical response of alu- 
minum under compression we calculated both the phonon 
dispersions of the fee phase and the cubic elastic con- 
stants (these were derived from the sound velocities, i.e. 
the slope of the phonons dispersions around T) as a func- 
tion of pressure. Experiments have shown that the fee 
phase may exist at pressures above the transition pres- 
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uncertainty^ At 222 GPa the predicted value of the c/a sure either as a super-pressurized phased or as a two 



ratio differs from experiment by 1%. Since experiments 
observe a region between 217 and 260 GPa in which the 
fee and hep phases coexist, this discrepancy is made more 
reasonable considering that the system might be out of 
equilibrium. We also show in Fig. [2Jd the equilibrium 
volumes for the different phases as a function of pres- 
sure. Volume differences between the phases are -0.055, 
-0.104, and -0.040 A 3 for the fee — > hep, fee — ► bee, and 
hep — > bee transitions respectively, corresponding to vol- 
ume changes of 0.6(6)%, 1.4(0)%, and 0.6(1)%. We note 
that in our phase sequence, and in those discussed in the 
literature^^ only the fee, hep, and bee phases are con- 
sidered. As a brief self-check we performed variable-cell 
relaxations at 1000 GPa using a four-atom unit-cell and 



phase region^ Therefore the mechanical properties of 
fee aluminum at pressures beyond the equilibrium tran- 
sition pressure, and any mechanical instabilities that may 
lead to mechanical failure are relevant to high pressure 
experiments. 

We calculated the phonon dispersions up to 1150 GPa 
using DFPT. For aluminum this method has been shown 
to accurately reproduce experimental values at P=0j^ 
Our calculations were performed with the 3e GGA pseu- 
dopotential, but compared at selected points in the Bril- 
louin zone with lie GGA calculations at pressures up to 
1200 GPa. As shown in figure [3] the discrepancies be- 
tween the 3e and lie results are at most of the order 
of 3-4% at the highest pressure, and much smaller be- 
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FIG. 5: Elastic constants as obtained from the sound veloci- 
ties along the [100], [110], and [111] directions. According to 
the Born criterion the fee phase becomes mechanically unsta- 
ble when Ac = cn — C12 = 0; this occurs around 725 GPa. 



low that. The phonon dispersions are shown in Fig. 0^ 
and we highlight the appearance, with increasing pres- 
sure, of a distinct softening of the lowest energy mode 
in the [110] direction. This starts to become evident at 
approximately 400 GPa, and is complete at 725 GPa, as 
highlighted in Fig. Since the slope of the dispersion 



curves is directly related to the elastic constants (Eq. [l| , 
we can extract the stiffness tensor from the vibrational 
modes near T; in our case the Born^ 8 - criterion for stabil- 
ity is 



1 dE 
H dki 



10 



Cn - C12 > 



(1) 



As Fig. O shows, the stiffness against shear deforma- 
tion, Ac = c\\ — C12, decreases above 400 GPa and goes 
to zero around 725 GPa, resulting in a Born 28 instability. 
These results complement existing studies of the proper- 
ties of bulk aluminu m 11 ! 13 and suggest another shear fail- 
ure mode, supporting previous studies suggesting shear 
failure modes to be inherent to bulk aluminum , 13 ' 15 More 
advanced treatments of mechanical stability including 
such effects of anharmonic modest and internal shear 
stresses created by loadin g 29 ' 30 ! 31 need to be considered 
in relation to the specific experimental setup before reli- 
able maximum stable pressures can be definitively deter- 
mined. 
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